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Abstract 

Highlights in the search for nonconventional (non qq) meson states are the 7Ti(1400) 
and 7Ti(1600) exotic candidates. Should they exist, mounting theoretical arguments 
suggest that they are tetraquark molecular resonances excitable by meson rescat- 
tering. We report a new tetraquark calculation within a model field theory ap- 
proximation to Quantum Chromodynamics in the Coulomb gauge supporting this 
■ conjecture. We also strengthen this claim by consistently contrasting results with 

exotic state predictions for hybrid (qqg) mesons within the same theoretical frame- 
work. Our findings confirm that molecular-like configurations involving two color 
singlets (a resonance, not a bound state) are clearly favored over hybrid or color- 
f — , ■ exotic tetraquark meson (qqqq atoms) formation. Finally, to assist needed further 

experimental searches we document a useful off-plane correlator for establishing the 
structure of these exotic systems along with similar, but anticipated much narrower, 
states that should exist in the charmonium and bottomonium spectra. 



The existence and understanding of exotic (non qq and qqq) hadrons is one of 
the few remaining closures to the standard model. Such states are expected 
according to Quantum Chromodynamics (QCD) and are of intense experimen- 
tal interest. In the light quark sector, there are two solid isovector candidates 
at 1.3-1.4 and 1.6 GeV, 7^(1400) and 7^(1600), each having J PC = 1~+ [IK] . 
Although doubts exist about the 1.4 GeV candidate [3J, a new analysis [I] 
supports it. The signature is observed as a p-wave resonance in the 77710 sys- 
tem and it therefore has odd parity P but even charge conjugation C yielding 
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the quantum numbers J PC = 1 h . Since a qq state with orbital L and spin 
S coupled to J = 1 must have P = (— 1) L+1 and C = (— 1) L+S , this state is 
clearly exotic. 



Assuming these states exist, the theoretical situation is even more controver- 
sial. The debate concerning their structure is among four possible scenarios: 
1) a hybrid (qqg) meson; 2) a tetraquark atom (qqqq involving intermediate 
color states that are not singlets); 3) a tetraquark molecular bound state of 
two conventional mesons; 4) a tetraquark molecular resonance (qqqq involving 
two intermediate color states that are singlets but not observed mesons). All 
four scenarios can produce J PC = 1~ + states but the first two are more exotic 
since the tetraquark molecule is color equivalent to a conventional meson- 
meson two-body state (see Fig. [I]). 



Lattice results, now performed with more realistic lower quark masses [3] , have 
focused upon the hybrid scenario but find that hybrid correlators can only 
produce a state as low as 1.9 to 2.1 GeV. A more recent lattice calculation 
[6] claims to find two hybrid meson masses below 2 GeV, however they use 
an ad-hoc extrapolation that is quadratic in the pion mass which increases 
uncertainties. Similarly, the lightest Flux Tube model 1 h predictions [TfgfD"] 
are also near 2 GeV spanning the region of 1.8 to 2.1 GeV. This is consistent 
with agreement among other model approaches, based on either the concept 
of constituent gluons [10] or field-theory calculations (see below) generating 
mass-gaps [TT], that the lightest hybrid mesons with just one constituent gluon 
should be somewhat heavier than the 1.4, 1.6 GeV experimental candidates. 
Finally, both well-established spin [12J and flavor |13] selection rules indicate 
that the above mentioned rj7r signature cannot be due to a hybrid meson decay. 
Therefore it would appear that the two tti states can not be theoretically 
explained as hybrid mesons. 
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Fig. 1. Four independent tetraquark color schemes. One is a singlet-singlet molecule 
while the other three are more exotic atoms (octet and two diquark schemes). 
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In search of other explanations, a potential model lifetime calculation p3] has 
ruled out a molecular ?7(1295)7r or ^(1440)7r bound state. However, it has been 
shown that meson-rescattering, specifically in the tjtt channel, could produce 
a resonance with this signature [T5JTB] - Also, it has been suggested [T7j that 
the 1.6 GeV resonance could be interfering with a background to produce the 
1.4 GeV structure. 

Summarizing the status of this situation, while the 7Ti(1600) can not be firmly 
precluded as a hybrid meson or exotic tetraquark, the 7Ti(1400) seems explain- 
able only as a molecular resonance excited by meson rescattering. The purpose 
of this work is to confirm the latter by a new theoretical analysis which also 
predicts that the 7r 1 (1600) is not a color exotic system. Our formalism, referred 
to as the Coulomb Gauge Model (CGM), has been successfully established in 
both the quark [T8lll9f20] and gluon [21] sectors and is based upon the exact 
QCD Hamiltonian in the Coulomb gauge [22J given by 



Hqcd = H q + H g + H qg + H c (1) 

H q = J dx^(x)[-ia ■ V + Pm]V(x) (2) 

H g =- jdx[j- l ll a {x) -JU a {x) + B a (x) • B a (x)] (3) 

H qg = g J dx J a (x) ■ A a (x) (4) 

Hc = ~y J dxdyp a (x)J- l K a \x, y)Jp\y) . (5) 



Here g is the QCD coupling, \1/ the quark field with current quark mass m, 
A a the gluon fields satisfying the transverse gauge condition, V • A a = 0, 
a = 1, 2, ...8, II a the conjugate fields and B a the non-abelian magnetic fields, 
B a = V x A a + \gf abc K b x A c . The color densities, p a (x) = ^(x)T a ^(x) + 
f abc A b (x) ■ H c (x), and quark color currents, J a = ^ (x)ctT ai &(x), entail the 
SU C (3) color matrices, T a = 4r, and structure constants, f abc . The Faddeev- 
Popov determinant, J = det(Ai), of the matrix M. = V ■ D with covariant 
derivative D ab = 5 ab V — gf abc A c , is a measure of the gauge manifold curvature 
and the kernel in Eq. (5) is given by K ab (x,y) = (x, a\A4^ 1 V 2 Ai^ 1 \y, b). 
The Coulomb gauge Hamiltonian is renormalizable, permits resolution of the 
Gribov problem, preserves rotational invariance, avoids spurious retardation 
corrections, aids identification of dominant, low energy potentials and does not 
introduce unphysical degrees of freedom (ghosts) [23J. To make the problem 
tractable, the Coulomb instantaneous kernel is approximated by its vacuum 
expectation value, yielding an effective potential field theory, Hqcd — > H^ CD 

H c - H$ = ~j dxdyp a (x)V(\x - y\)p a {y) , (6) 

with confinement described by a Cornell potential, V(r) = — — + or, where 
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the string tension, a = 0.135 GeV 2 , and a s = 0.4 have been independently 
determined from conventional meson studies [T8|19|20j within the same field 
theory approach. We also use the lowest order, unit value, for the Faddeev- 
Popov determinant in the H g term and treat the H qg interaction using per- 
turbation theory. Lattice data confirms the Cornell potential form between 
static sources [21] and further provides the scale of the gluon mass gap [25] . 
that is needed to fit a counterterm in the gluon gap equation. The remaining 
parameters are the reasonably well known current quark masses at some high 
energy scale where the mass function runs perturbatively. The quark sector 
Hamiltonian then takes the form of the Cornell coupled- channel model |26j . 
We note that three-body forces [2~Tf28] are omitted, however based upon suc- 
cessful three-body applications [21] we submit the CGM should capture the 
dominant features of a multi-parton spectrum. 

Before presenting our tetraquark results we highlight our recent hybrid me- 
son calculation [29|. Since the gluon carries a color octet charge, the quark 
and antiquark are also in a color octet wavefunction with elements (they 
repel each other at short distance). In the hybrid rest frame there are 
two independent three-momenta q + = S^ 2 , q_ = q — q and one dependent 
g = —q — q = — 2q + . The leading hybrid Fock space wavefunction can there- 
fore be constructed from the respective quark, anti-quark and gluon quasipar- 
ticle operators B{ iCi (q) , D{ 2C2 (q) and afi(g) 

An angular momentum expansion for the lightest 1 h state reveals a required 
p-wave excitation in one of the two orbital wave functions. Consult Ref. [21] 
for complete details. Significantly, in agreement with earlier findings [11], the 
predicted hybrid masses are about 2 GeV for the ground state (parity +) 
quadruplet and 2.2 and 2.4 GeV, respectively, for the first 1~ + exotics. The 
repulsive nature of the short range qq interaction potential kernel and the large 
gluon mass gap are responsible for these large masses. We also performed a 
parameter sensitivity and error analysis study and concluded there was no 
model possibility to lower one of these states near the 1.6 GeV candidate and 
therefore rule out this, and even more clearly the 1.4 GeV, state as a hybrid. 

Returning to our thrust, we report results for qqqq spectroscopy (see also a 
preliminary study [30]). This system was first investigated in the bag model 
[5T] and then more extensively by Ref. [32] with subsequent potential model 
applications reported by Refs. [3"3"f34"] . While these studies have some similarity 
to our model, we submit our results are more robust since our approach is 
much more comprehensive, has many QCD elements with no new parameters 
to be determined and employs a realistic potential kernel extracted from lattice 
gauge theory. The leading tetraquark Fock space wavefunction is [3"51I3"6] 
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^jp C) _ ff fdq A dq B d Ql JPC , , 

^g<ft(9i)^U(9 a ) s U(9s)^U(94)l n >- ( 8 ) 

In the cm there are three independent momenta which we take to be 
q A = 91 ~ 92 , q B = <73 2 <?4 and qj = q ^ qA — SilSb. (g^ g 3 f or the quarks and q 2 , 
q 4 for the anti-quarks). The color matrices RqIcI yielding color-singlet wave- 
functions follow from the SU C (3) algebra depicted in Fig. [TJ The radial part 

2 2 2 

of H?£)*)«(Qa, <?b, <?/) is chosen to be a gaussian, exp(-^ - J| - ||), with 
variational parameters a^, and aj for s-wave states and a gaussian multi- 
plied by qf/af (i = A,B,I) corresponding to orbital Li = 1, when treating 
p-wave states. Note, as for hybrid mesons, the ground state tetraquark multi- 
plet has positive parity and that constructing 1~ + exotics requires one of the 
three orbitals to be a p-wave. Using the variational principle, the tetraquark 
mass is then given by 



Mjpc < 



l^,JPC 
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n QCD 
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M self + M qq + M m + M qq + M, 



annih 



(9) 



Contributions to the Hamiltonian expectation value are summarized in Fig. [2] 
and correspond to 4 self-energy, 6 scattering, 4 annihilation and 70 exchange 
terms, each of which can be reduced to 12 dimensional integrals that are eval- 
uated in momentum space. Because of the computationally intensive nature of 
this analysis, the hyperfine interaction was not included. Complete expressions 
will be given in another publication, but as an example note the annihilation 
contribution, not possible in standard quark models, is 
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Fig. 2. Equal-time diagrams for the expectation value of the model Hamiltonian. 
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12 3 V(q x + q 2 )u\, ( 9l + fc)iv (g 2 - fe) x 



(27T) 

^(^MgO^fw^i' 9 2 > ^^^fai + fc> <l2 - k, q 3 ), (10) 



involving Dirac spinors u\ 1 and v\ 2 . This raises the mass of the isoscalar states 
relative to the states with higher isospin, which would otherwise typically be 
heavier due to the exclusion principle applied to equal-flavor quarks. 

Performing large-scale Monte Carlo calculations (typically 50 million samples), 
has conclusively determined that the molecular representation (i.e. singlet- 
singlet) produces the lightest mass for a given J PC . This is due to suppression 
of certain interactions in this color scheme from vanishing color factors for 
every parton pair which does not occur in the other representations. Also, 
there are additional, repulsive forces in the more exotic color schemes. Using 
fn u = Tn>d — 5 MeV, the predicted tetraquark ground state is the non-exotic 
vector 1 ++ state in the molecular representation with mass around 1.2 GeV. 
Figure [3] depicts the predicted tetraquark spectra for states having conven- 
tional and exotic quantum numbers in both singlet and octet color represen- 
tations. The quark annihilation interactions (qq — > g —>■ qq) in the I g q = 
channel generate isospin splitting contributions, up to several hundred MeV, 
in the octet but not singlet scheme as slightly illustrated in the figure. Isospin 
splitting is a consequence of a more proper field theory treatment, not present 
in conventional quark models. The annihilation interaction terms are repul- 
sive, yielding octet states with 1 = 2 lower than the 1=1 which are lower 
than the 1 = 0. This is intuitively contrary to expectations that 1= 2 states 
are higher based upon the Pauli principle that identical quarks repel. The 
molecular states are all isospin degenerate and the lightest exotic molecule is 
an intriguing with mass 1.35 GeV that could be detected in a sophisti- 
cated p-wave analysis of an un spectrum. Because of the isospin degeneracy, 
there will be several molecular tetraquark states with the same J PC in the 
1 to 2 GeV region. Further, these states can be observed in different elec- 
tric charge channels (different I z ) at about the same energy, which is a useful 
experimental signature. The lightest 1 h is predicted near 1.4 GeV which is 
close to the observed tt(1400), suggesting this state has a molecular resonance 
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Fig. 3. Tetraquark 1=1 singlet (molecule) and octet (atom) schemes spectra. 
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Table 1 

Selected tetraquark molecular exotic states. For qq pairs, Si and Lj, i = A,B, are the 
total spin and orbital angular momentum and Lj is the orbital angular momentum 
between the pairs. The variables not shown are assumed to be 0. Units are GeV. 
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structure. The computed mass for 1 h states with more exotic octet color 
configurations are all above 2 GeV. This is consistent with model predictions 
[29] for exotic hybrid meson (qqg) 1 _+ states also lying above 2 GeV due to 
repulsive color octet quark interactions. Finally, for any J PC state, including 
the 1~ + , the computed masses (not shown) in both the triplet and the sextet 
diquark color representations are all heavier than in the singlet representation 
and comparable to the octet scheme results. Our predictions for the lightest 
exotic molecules are given in Table [TJ 

It is interesting to document the current quark mass dependence of our results. 
This is illustrated in Fig. HI where the calculated rest mass corresponding to 

Tetraquark and two meson masses 
Quark mass dependence 

, + — 1 4000 

: 

.* - 3000 
// 
•/ 
'/ 
/' 

// : 
A // - 2000 

-•'/ 



- 1000 



u i io ioo ' iooo 

m q (MeV) 

Fig. 4. Dotted lines, from bottom to top, are the 7T7t, ttt], ttt]' ,nr] c observed thresholds, 
respectively. The two other plots are the pseudoscalar m^+Mgg model predictions in 
the (chiral respecting) RPA and (chiral violating) TDA. The squares and diamonds 
correspond respectively to predicted scalar isoscalar tetraquark masses in the color 
singlet-singlet and triplet-triplet schemes. Note the two-meson calculation [19] has 
a = 0.18 Gev 2 , a s = 0, whereas our results use a = 0.135 Gev 2 , a s = 0.4, but the 
difference among both sets is known to be small for the ground state. 
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Table 2 

Predicted masses of hidden-charm exotic 1 h mesons (molecular configuration). 



m c = 1.2 GeV 


m g 


M mq - 


ccqq 





4.04 GeV 


ccqq 


110 MeV 


4.10 GeV 


ccqq 


150 MeV 


4.15 GeV 



different four quark systems (tetraquarks and two mesons) is displayed with 
one qq pair being held fixed at a low 1-5 MeV while the other pair varies from 
1 to 1300 MeV. Note the predicted ccuu tetraquark mass is relatively closer 
to the two meson decay threshold indicating a narrow decay width to the irr] c 
and 7]f] c channels (more so for bb, not displayed). Also note a quark-meson 
exchange model [37] claims some J = 1 partners will be absolutely stable. 
The lowest-lying hidden-charm exotics are predicted in Table [2] and can be 
searched for in hidden charm decays (p-wave r) c 7To, r\ c r\ for example). 

Figure [5] displays the sensitivity of the scalar tetraquark mass to the variational 
parameters for different color schemes. The color singlet- singlet configuration, 
at 1.28 GeV, is clearly the lightest, in agreement with previous predictions 
[38] . Similar results are obtained for other J PC configurations. This suggests 
that the lightest tetraquark states, including the 7Ti(1400) should it exist, are 
best interpreted as molecules of color singlets (the "extraordinary hadrons" 
|39j ) and directly excitable via meson rescattering. 
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Fig. 5. Tetraquark mass of the ground state isoscalar ++ as a function of the 
intercluster variational parameter, for various color configurations. Note the more 
exotic color configurations are above 1.7 GeV. 
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Fig. 6. Hybrid (left) and tetraquark (right) meson probability density contour plots. 
The left horizontal axis is half the momentum of the constituent gluon, q + = g/2. 
The left vertical axis is the angle between the qq relative momentum q_ and the 
hybrid total angular momenta J. The right horizontal axis corresponds to the in- 
tracluster momentum qj and the angle between this momentum and J for the 
tetraquark. Angles in radians, momenta in GeV. 

For further dynamical insight, Fig. [6] details contour plots of the probability 
densities using the hybrid (left) and tetraquark (right) variational wavefunc- 
tions. Note the depletion of the wavefunction at low-momentum, reflecting the 
rising confining potential at large distance, and the significantly different par- 
ton momentum distributions between hybrids and tetraquark systems. This 
difference in momentum distributions will produce distinct decay signatures 
which can be utilized to identify hybrids and tetraquarks as we now discuss. 

It has been proposed jlO] that 1 h exotic hybrid mesons decay preferentially 
to a meson pair in a relative s-wave, where one of the mesons is a p-wave 
(axial) meson. However that prediction was based on the Flux Tube model 
and not known to hold exactly in any limit. We have recently [29J detailed a 
less model dependent decay signature based upon the Franck-Condon principle 
of molecular physics, which predicts that the momentum distribution of decay 
products parallels the internal momentum distribution of the parent meson. 
This is an exact statement in the infinitely heavy quark limit, however it is 
also useful for physical quarks. Our previous application [29] was designed 
to distinguish qq from qqg mesons and is now generalized to the four-body 
problem. For two and three-body systems, in the cm partons' momenta are 
restricted to a plane, however for tetraquarks there are off-plane degrees of 
freedom which can be exploited as an identification signature (see the cartoon 
in Fig. [7]). One observable sensitive to this is an exclusive measurement of 
decay to four mesons. By then determining their momenta p i=1 2 3 4 in the 
cm of the exotic candidate, kinematic cuts can be applied for every two and 
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Fig. 7. Differences between internal in-plane three-body hybrid (top left) and off- 
plane tetraquark (top right) meson decays. From the Frank-Condon principle the 
internal momentum distribution is reflected in the momentum distribution of the 
final products yielding an off-plane final state momentum observable, II, that is very 
different for qqg (bottom left) and qqqq (bottom right) mesons. 

three meson groups to eliminate intermediate two and three-body resonances 
that could confuse the analysis. From the pure (and much reduced) four-body 
decay sample count, one can construct the off-plane correlator, II, for any three 
momenta which is the volume of the parallelepiped they form in momentum 

space, 



n(pi,P2,p 3 ,p 4 ) 



((Pi XP 2 )-P 3 ) 2 



Pi x P2IIP2 x P3IIP1 x P 3 ||Pi x p 4 ||p 2 x p 4 ||p 3 x p 4 \ 



which has been normalized by the area of the six faces. This dimensionless 
off-plane correlator is useful for the very different scales involved in light and 
heavy quark physics, positive definite and invariant under permutation of the 
four momenta. In the cm system only three momenta are linearly independent. 
Taking them equal and along the edges of a cube yields II = 8 -1 / 4 = 0.59, 
with a more general maximum value close to 0.707. Therefore the value of the 
correlator event by event will be a random variable distributed between and 
0.707 and, with sufficient statistics, one could deduce the internal structure 
of the decaying resonance (bottom right plot of Fig. [7]). This could then be 
compared to four-meson decays of established conventional and hybrid meson 
benchmarks, as distinguished using the procedure outlined in Ref . , which 
would be radically different (their II volumes collapse to in the heavy quark 
limit). Finally, we are currently calculating tetraquark decay widths to com- 
pare with both conventional and hybrid meson decays to two and multi-meson 
final states. These will also reflect hadronic structure differences and therefore 
aid exotic identification. Results will be reported in a future communication. 
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In conclusion, we submit that the observed 1 _+ exotics below 2 GeV are 
not color exotic hadrons but rather somewhat more conventional tetraquark 
molecular resonances involving color singlet qq pairs. Our results also agree 
with lattice and other models predicting that color exotics (hybrid mesons, 
octet-octet, sextet-sextet and triplet-triplet tetraquarks) will have masses near 
and above 2 GeV. Finally, we have utilized the Frank Condon principle and 
the distinctive quark momentum distribution in a tetraquark to propose a new 
off-plane correlator measurement for identifying exotic hadrons. 
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